Fitting Formula for Flux Scintillation of Compact Radio Sources 



J. Goodman 
Princeton University, Princeton, NJ 08544 
j eremyOastro . princeton . edu 
and 
R. Narayan 

Harvard- Smithsonian Center for Astrophysics, Cambridge, MA 02138 
narayanOcf a . harvard . edu 

ABSTRACT 

We present a fitting function to describe tlie statistics of flux modulations 
caused by interstellar scintillation. The function models a very general quantity: 
the cross-correlation of the flux observed from a compact radio source of finite 
angular size observed at two frequencies and at two positions or times. The 
formula will be useful for fitting data from sources such as intra-day variables 
and gamma-ray burst afterglows. These sources are often observed at relatively 
high frequencies (several gigahertz) where interstellar scattering is neither very 
strong nor very weak, so that asymptotic formulae are inapplicable. 

Subject headings: scattering - radio continuum: general - methods: data analysis 



1. Introduction 

Ever since its initial discovery in the signals received from radio pulsars, interstellar 
scattering and scintillation has proved to be a useful tool in radio astronomy (see Rickett 
1990, 2001; Narayan 1993; Hewish 1993 for reviews). It provides unique information on 
small-scale turbulent fluctuations in the ionized interstellar medium (e.g., Armstrong, Cordes 
& Rickett 1981; Wilkinson, Narayan & Spencer 1994) and on the sizes of compact radio 
sources. The latter application has contributed to progress in several areas of astrophysics: 
radio pulsars (e.g., Roberts & Abies 1982; Cordes, Weisberg & Boriakoff 1983; Wolszczan & 
Cordes 1987; Gwinn et al. 1997), intra-day variables (e.g., Rickett et al. 1995, 2001; Lovell 
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et al. 2003) and gamma-ray burst afterglows (e.g., Goodman 1997; Frail et al. 1997; Taylor 
et al. 1997). 

The fundamental quantity in scintillation theory is the correlation of flux variations of a 
distant source observed at two positions separated by transverse distance r on the observer 
plane. Since each flux is the square of the local electric fleld, the quantity of interest is a 
fourth order correlation of the electric field (see eq. 2). To interpret scintillation data, wc 
need to calculate the expectation value of this fourth moment as a function of the strength of 
the turbulent fluctuations in the scattering medium, the shape and size of the radio source, 
the observing frequency, etc. 

The scattering of radio waves in the interstellar medium may be characterized by two 
dimensionless parameters (eq. 6): U which measures the strength of the scattering, and a 
which describes the power spectrum of the fluctuations. The turbulence in the interstellar 
medium^ appears to be well-described by a Kolmogorov scaling, which corresponds to a = 
5/3. Hence, the scattering medium along any line of sight is determined by just U ; however, 
U varies from one direction to another, and also with frequency in a given direction (see eq. 
19). Analytical results are available for the flux correlation function in the hmit of both very 
weak scattering (U <S 1) and very strong scattering {U » 1) (e.g., Goodman & Narayan, 
1995; Rickett 1990; and references therein). In the latter regime, the scintillation is known 
to occur on two very different scales, determined by diffractive and refractive effects. 

While a great deal of interesting work on interstellar scintillation has been done using 
asymptotic results, many observations correspond to the difficult intermediate regime where 
U is neither -C 1 nor 3> 1; no analytical resiflts are available in this transition regime. For 
typical high lines of sight through the interstellar medium at high Galactic latitudes, the 
transition regime occurs at radio frequencies ~ 5 — 10 GHz, a frequency of much interest 
for both intra-day variables and gamma-ray burst afterglows (but less so for pulsars which, 
because of their steep radio spectra, are usually observed at lower frequencies where asymp- 
totic strong scattering results apply) . Some very approximate formulae have been proposed 
for the transition regime (e.g.. Walker 1998). However, for accurate results, one has to re- 
sort to numerical computations of the scintillation correlation function, which is technically 
challenging and has been rarely attempted. 

We describe in this paper numerical computations of the scintillation correlation function 
for a wide range of values of the scattering parameter U, spatial separation r, source size 
(deflned in eqs. 15-18) and frequency difference 77 (deflned in eq. 5). Using the numerical 



^We use the term "turbulence" as a short-hand for both true dynamic turbulence as well as passive 
density irregularities that may not exhibit turbulent dynamics 
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results we have developed a fitting function for the flux correlation that is vahd for all values 
of U. The function asymptotes to the appropriate analytical results in the limits [/ -C 1 
(very weak scattering) and U ^ 1 (very strong scattering) and agrees well with numerical 
results in the transition regime in between. 

In §2 we introduce the basic fourth order moment of interest to flux scintillation obser- 
vations, and in §3 we explain how we numerically compute this quantity. In §4 we present 
our fitting formula, which agrees well with the numerical results. We conclude in §5 with 
a brief summary. We present in Appendices A-C technical details, including some asymp- 
totic results useful as checks of our fitting functions, and a discussion of parabolic arcs in 
secondary spectra. 



2. Thin-screen theory 

In the thin-screen approximation, one imagines that the turbulence affecting radio-wave 
propagation is concentrated in a narrow layer perpendicular to the line of sight. This is 
sometimes not far from the truth, although it would usually be more accurate to assume a 
number of scattering screens; but the opposite limit of homogeneously distributed turbulence 
is probably the exception rather than the rule. The observer lies in a parallel plane at distance 
^screen from the scrccu, and the source is behind the screen at distance -^source from it. Waves 
propagate freely from the source to the screen, where they suffer a phase shift 0(x) that 
depends upon the two-dimensional position x within the screen, and thence travel freely 
again to the observer's plane. The distortion of the phase fronts at the screen leads to 
modulations of the flux on the observer's plane, through a combination of refractive focusing 
and diffractive interference. It is useful to define an effective distance 

\ ■^screen ■^source / 

from observer to screen. This allows us to analyze the problem as if the source were infinitely 
distant and the wavefronts were planar before encountering the screen. 

We follow the notation of Goodman & Narayan (1989, henceforth GN89), who discuss 
the mutual coherence function 

r4(b; r; ly,, v^) = Ie Qb, E* (-\, E (v + ^b, E* (r - ^b, ^^2) ) , (2) 



where E{r, u) is the electric field due to a source of unit strength measured at vector position 
r on the observer's plane (a plane perpendicular to the line of sight to the source) at radio 
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frequency v. As written above, r4 is useful for studying the correlation between measure- 
ments of visibility made with two two-element interferometers having the same baseline b 
but phase centers separated by r, and operating at different frequencies z/i and z/2. In this 
paper, we are interested in flux correlations rather than visibilities, so we set b — > 0. Posi- 
tions and r then represent two positions at which the flux is measured. The measurements 
may be done by two telescopes simultaneously, or by the same telescope at different times. 
In the latter case, the interstellar turbulence responsible for phase changes in E is assumed 
"frozen" on timescales of interest, so that r = Vj_i, where Vj_ is the effective transverse 
velocity of the hue of sight through the interstellar medium. The angle brackets denote an 
ensemble average over realizations of the turbulence. 

Appendix A derives formula (A4) for the flux correlation r4(0; r; z/i, U2) in terms of the 
phase structure function 

L>(Ax) = ([0(Ax,A)-0(O,A)n 
evaluated at the geometric mean wavelength 

A = Va^ (3) 

of the two observing wavelengths Ai,2 = c/ ^1,2-1 and in terms of a Presnel scale deflned with 
the arithmetic mean, 

{\1 + \2)Z^"' 



and a dimensionless frequency difference 



47r 



(4) 



V^^^-^ : 0<|r7|<l. (5) 

In this paper, we restrict our attention to power-law structure functions 

^(Ax) = C/(^y, (6) 

and unless otherwise speciflcd, wc assume a = 5/3 (the "Kolmogorov" value). Thus C/ is a 
dimensionless number describing the strength of scattering; it depends upon the observing 
wavelength (see eq. 19) as well as the intrinsic amplitude of the turbulent electron-density 
power spectrum, since D oc A^ as a consequence of the plasma dispersion relation. Henceforth 
we choose units of length such that rp = 1. Also, we normalize the mean flux to unity at all 
frequencies, (-F(r, i/)) = {\E{r, z/)p) = 1. Writing W rather than r4 for the flux correlation, 
we have 

W{r,ri) = U{0■,r■,u,,U2) = {\E{0,u^)\'\E{r,u,)\')^ y exp(ir • q)iy(q, 77), (7) 



- 5 - 



W{q,rj) — exp 



ds exp(iq ■ s) exp 



F{q,r)) = U 



(1 - r?)"+^ + (1 + r?)"+^ 



G{s, q,rj) = U [\s - ?7q|" + |s + r)q\" - |s - q|" - |s + qp] . 

The quantity W{q^,r]) will be referred to as the "cross spectrum." For rj 
the two-dimensional spatial power spectrum of the flux, H^(q). 



(9) 
(10) 

it reduces to 



We know no closed-form results for the integral (8) except at uninteresting values of 

a, so it is necessary to evaluate the integral numerically in general. However, approximate 
analytic expressions can be obtained in the limits -C 1 (weak scattering) and [/ S> 1 (strong 
scattering), and these are given in Appendix B. They are useful to check the accuracy of the 
numerical evaluation of the integral (7), and also to guide the choice of functional form for 
the fitting functions presented in §4. 



3. Numerical methods 

The integral (8) for the flux cross spectrum Vr(q, 1]) proved to be challenging to estimate 
accurately over the full ranges of U, q, and 7] required. It is two-dimensional and extends to 
infinity. Worse, it is only conditionally convergent since G'(s, q, t]) ^ as s/q ^ oo. And 
the regions of the s plane that dominate the integral vary with the scattering regime. 

We experimented with several approaches before settling on the following. The inde- 
pendent variable s is resolved into its components parallel and perpendicular to q: sy & s^. 
The inner integration is taken over sy and is calculated by a version of steepest descent in 
the complex plane s^^ = z = x + iy. That is, at each value of s± required by the quadrature 
scheme for the outer integral, our code estimates 

j exp[iqz + G{z, s±,q,r])] dz, (11) 

c 

where C is a contour which in the first instance is the real axis [y = 0, x E {—oo, oo)] but 
which is moved into the upper half plane, where the first exponential in (11) vanishes as 
y — > oo since q = |q| > 0. The contour must not, however, be dragged across any of the 
singularities of 

G{z,s^,q,v) = u{[{z-vq? + slY^'+[{z + vqr + slY^' 

-[iz-qr + slf'+[{z + qr + sir^'}, (12) 
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which has algebraic branch points at y = s_l and x e {±?, ±^?}- A branch cut must be 
drawn from each point. Unless Aa is integral, these cuts must extend to infinity rather than 
join the branch points to one another. We chose to draw the cuts upward parallel to the 
imaginary axis. This is convenient because [{z — of + s^]"/^ is evaluated as 

exp{(Q;/2)Log [{z - of + s\]] , 

and complex library functions of scientific programming languages usually put the branch 
cut of Log on the negative real axis, which corresponds to the cuts in the z plane that we 
have described. For pedagogical reasons, we wrote the code in C, although that language 
has only recently supported complex arithmetic. Recent versions of the GNU gcc compiler 
and mathematical library functions performed very reliably and were almost as convenient 
to use as their Fortran counterparts. 

A typical contour of integration is shown in Figure 1. Integration starts from x ~ 0, 
since the complete integral (11) is twice its righthand half, but at y > 0, since there is no 
singularity to pin the contour to the origin. The contour is drawn adaptively, differentiating 
the argument of the exponential at each step to determine the direction in which its real 
part becomes negative most rapidly and its imaginary part is constant, until a branch cut is 
encountered. The contour then follows the lefthand side of the cut downward to the branch 
point, circumnavigates it, and then continues by steepest descent. Special care must be 
taken when q is small since there can be near-cancellation between the contributions from 
the cuts eii X — rjq and x — q; this is handled by applying a Romberg quadrature scheme to 
the sum of the integrands a.t z = {rjq — 0"*", y) and {q — 0"*", y). Care must also be taken when 
|s| 3> g because the terms of G nearly cancel; this is handled by expanding G in powers of 
q/z. 

The numerical estimate of (11) becomes the integrand for integration over s±. This 
integrand decreases exponentially at large s± because the branch points of the inner integral 
move far from the real axis where their contribution is suppressed by ~ exp(— gs^). We 
use a straightforward midpoint quadrature along the real axis in the auxiliary variable t = 
sinh~^(sj^/so), where Sq is a scale chosen according to the values of q and U. As is well 
known, numerical quadratures of smooth functions over infinite or periodic domains with 
schemes that use uniform steps and weights, such as midpoint or Simpson's rule, converge 
faster than any power of stepsize. So there is no point in using a higher-order method for 
the outer integral. 

The present method for evaluating W{q) is efficient for moderate to large values of U 
and q because the s\\ integrand decreases rapidly along the steepest-descent path, but it has 
difficulty with very small Uq°'. Fortunately, the analytic weak-scattering approximation is 
then very accurate, so the code uses that approximation when Uq°' < 10~^. One test of 
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Fig. 1. — Integration contour (solid line) and branch cuts (dashed) for a = 5/3, C/ = 4, 
rj = 0.2, q^l, s± = 0.1933. 
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the code is that the numerical results match smoothly onto the analytic ones in the weak- 
scattering regime. We have tested the code also against the asymptotic results of the previous 
section for strong scattering, t/ 3> 1. And we have tested it against an old code valid only 
for 7] = that uses a completely different strategy (based on integrating s in polar rather 
than cartesian coordinates; see Goodman & Narayan 1985). 

After tabulating VF(q; U, rj) for given values of U and rj, we obtain the flux correlation 
function W{r) by numerical calculation of the Hankel transform 

oo 

W{U, r,v)-^j Mqr)W{q; U, rj) qdq , (13) 



where W(r) and W{q) depend on the absolute values of r and q only. (This would not be true 
of the spectrum of phase fluctuations on the scattering screen were anisotropic.) Equation 
(13) is for a point source. In the general case, when the source observed at position 
has a circularly symmetric normalized intensity profile Si{r) with root mean square size r^i 
and that observed at position r has profile S2{r) with size rs2, the flux correlation function 
becomes 

oo 

W{U,r,rsi,r,2,v) = ^J Mqr)S,{q)S2{q)W{q;U,rj) qdq , (14) 



where Si{q) and ^2(5) are the Hankel transforms of the source profiles: 

00 

SiM = J Joiqr)Si,2ir)rdr. 


In this paper, we consider two different source profiles, gaussian and tophat. The normalized 
profiles are 

Gaussian : S(r) ^ ^e-^'^''^^' ; Siq) ^ e-^'^'^'^'^' (15) 

Tophat : S{r) = ^H{rs-T); S{q) = —J,{qrs) , (16) 
7rr| qrs 

where is a measure of the source size, and H{x) = 1 if x > 0, H{x) = if x < (Heaviside 
function) . Note that the root mean square radius of the source Tj-ms is related to the nominal 
size rg as follows: 

rrms = rs (Gaussian), r^^s = rs/\/2 (Tophat). (17) 
The ratio Trms/'^s has an influence on the constants Ci and C2 defined in §4.2. 
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4. Fitting Function 

In developing a fitting function for the correlation W{U,r,rsi,rs2,i]) equation (14), 
we have found that it is not necessary to consider the dependence of W on the two sources 
sizes Tgi and independently. Rather it is sufiicient to consider a single effective size 

rs^{rl, + rl,y'\ (18) 

Thus, the most general scintillation flux correlation function that we consider in this paper 
is l^(t/, r, Tg, r]), where U measures the strength of the phase fluctuations at the scattering 
screen at the geometric mean wavelength A (eqs. 6 and 3), r is the distance in Presnel units 
Tp (eq. 4) between the two points at which the fluxes are measured, is the effective source 
size in Presnel units as projected at effective distance z (eq. 1), and 77 is the dimensionless 
frequency difference between the two observations (eq. 5). All the results presented here 
are for a Kolmogorov spectrum of fluctuations in the scattering screen (a = 5/3). For this 
choice of a, the parameter U varies with wavelength A as 

/ \ \ (4+a)/2 

U(X) ^ (3^) , (19) 

where Aq is the wavelength at which U — 1. Thus, the scattering power of the screen is 
uniquely deflned by the wavelength Aq. 

Using the code described in the previous section, we have calculated numerical values of 
W(U, r, Tj) over a wide range of values of U, r, and rj, and we have developed a fltting 
function that approximates the numerical results. 



4.1. Mean Squeire Flux Variations for a Point Source 

In the limit when r = = 77 = 0, the quantity 

Wo{U) = W(U, r = 0, r, = 0, 77 = 0) (20) 

measures the mean square flux variations due to scintillation as a function of the scattering 
strength U. The panel on the left in Figure 2 shows a numerical evaluation of this function. 
An accurate fltting function is given by 

0.7729[/ 1 + 0.4760[/~°-^ 

^ ~ 1 + 0.286t/ + 0.0860[/2 + 0.0550[/3 ^ 1 - IMU-^ + lO.lt/-2 ' ^ ' 

which has an error much less than 1% over all values of U (Fig. 2, right). The function 
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Fig. 2. — Left: The solid line is the correlation Wo{U) as calculated numerically, and the 
dashed line is the fitting function (hardly seen). Right: The difference between the calculated 
correlation and the fitting function. 
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Wo{U) crosses unity at 

Uo = 1.5874. (22) 

We use this value of U to distinguish between the two major regimes of scattering. When 
U < Uo, we consider that we are in the weak scattering regime and seek to explain the 
scintillation properties using a single scale, the weak scale. For U > Uo, on the other 
hand, we consider that we are in the strong scattering regime and allow for two scales, one 
for diffractive scintillation and one for refractive scintillation. These are explained in the 
subsections below. 

4.2. Weak Scattering {U < Uo) 

In this regime, we write the fitting function for the general correlation W{U, r, Ts, rj) in 
the form 

W{U, r, r„ r,) = Wo{U)FrF,F^, (23) 

where 



Fr 


- exp [-(r/i?i)5/3] , 


(24) 




= [l + {rjR,y-'Y'^ 


(25) 




= [1 + a,{U)C'/' + a^mey', 


(26) 


c 


= 77/(1-77), 


(27) 


Ri 


= i?3 + 3.5477^-^7^3 


(28) 


R2 


= i?4 + 3.5477l•lVi^l•°^ 


(29) 


Rs 


= {Rl + C,rlY^\ 


(30) 


R4 


— C2R5, 


(31) 


R5 


= 1.15 - 0.260([//[/o)°'^^°, (weak scattering only). 


(32) 


ai(U) 


= [l.0l2 + (0.747C/)2]'/', 


(33) 


a2{U) 


= "(ll.9t/i-02)V (5.37^72-4) ']'^'. 


(34) 



The three factors Fr, Fg and describe the variation of the fiux correlation as a function 
of separation r, source size and wavelength difference 77. Each factor is defined such that 



it goes to unity in the limit when r = = 77 = 0, so that equation (23) reduces to equations 
(20) and (21) in this hmit. 

There are two constants, Ci and C2, in the above formulae. The values of these constants 
depend on the shape of the source. For the gaussian and tophat source models defined in 
§ 3, we find 

Gaussian : Ci = 1.37, C2 = 0.965, (35) 
Tophat : Ci = 0.829, C2 = 1-32. (36) 

All the other numerical constants in the fitting function are independent of source shape. 
We suspect that the above two source models are sufficient for most applications — the 
gaussian will serve for the majority of sources and the tophat is suitable for gamma-ray 
burst afterfiows (Sari 1998). In case one needs to consider other source shapes, we note that 
Ci and C2 scale approximately as 

Ci«.1.4(^y, C2^0.95('^), 

where rrms is the root mean square size of the image. Equation (17) gives r\ms/^s for the 
gaussian and tophat models. The scaling (37) is rather approximate because the best-fit 
values of Ci and C2 depend on more than just the second moment of the source. 



(37) 



4.3. Strong Scattering {U > Uq) 

When U > Uq wc are in the regime of strong scattering and the flux variations have 
contributions from both diffractive and refractive scintillation. In the limit when r = = 
77 = 0, we assume that the mean square flux variations Wo{U) consist of contributions 
Wd{U) and Wr{U), respectively, from diffractive and refractive scintillation. We take these 
contributions to be given by 

Wa{U) = 

The two quantities add up to give Wo{U), as they should. They are also so defined that 
when U ^ Uq, the refractive term vanishes, i.e., at the boundary between weak and strong 
scattering, we have only diffractive scintillation. Thus, the fitting function will be continuous 
across Uq if we match the weak scattering function (23) of the previous subsection with the 
diffractive term in the function (39) below. 



Wq{U) 



Wr{U) 



Wq{U) - 1 



(38) 
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With this motivation, we write the fitting function in the strong scattering regime as 

W{U, r, r„ r,) = Wd{U)FrFsF, + (39) 

where the first term on the right is the diffr active term and the second is the refractive term. 
We take the factors Fr, Fg and in the diffractive term to have the same forms as in § 4.2, 
with the sole exception that the scale is now given by 

i?5 = 0.964C/-°-^ + (0.890 - 0.964C/-°-^) (C/o/C/)^■^^ (strong scattering only). (40) 

The two expressions (32) and (40) for are equal ai U = Uq. Therefore, we have a perfect 
match at f/ = Uq between the fitting functions for weak scattering and strong diffractive 
scattering. 

The second term on the right in equation (39) describes refractive scintillation. Here 
we set 



TP 

± r,r 


= exp [-(r/it:6)'/^] , 


(41) 


^ s,r 


= [l + {rjR,y"]-\ 


(42) 


F 

^ V,r 


= [l + a,{U)C''^' + a,{U)CY\ 


(43) 


Rq 


= i?8+l5.v■^^ 


(44) 


Rv 


= i?9+l5.lv■^^ 


(45) 


Rs 


= {Ri,+cy:)"\ 


(46) 


Rg 


= C2R10, 


(47) 


RlO 


= 1.20 - 2.00t/°-3 + 1.73t/° ^ 


(48) 




= a,{Uo){Uo/Uf-''', 


(49) 




= 8.38 + [a2{Uo) - 8.38] {Uo/Uf^^^. 


(50) 



The constants Ci and C2 have the same values as before. 

4.4. Comparison with Numerical Results 

Figures 3, 4, 5 show the dependence of the flux correlation function W{U,r,rs,r]) on 
each of the three variables r, and 77, while keeping the other two variables fixed at zero. 
Both the exact numerical results and the fitting function are shown. 



-14- 





Fig. 3. — Solid lines show the correlation as a function of separation r for zero source size 
{vg — 0) and zero wavelength difference (r/ = 0) as calculated numerically. Dashed hues 
show the corresponding fitting function. Top Left: log(C/) = —0.25, — 0.5, — 0.75, — 1.0 
(from above). Top Right: log{U) = 0.5, 0.25, (from above). Bottom Left: log(C/) = 
0.75, 1.0, 1.25, 1.5 (from above). Bottom Right: log(C/) = 1.75, 2.0, 2.25, 2.5, 2.75, 3.0 
(from above). 
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Fig. 4. — Solid lines show the correlation as a function of source size for zero separation 
(r = 0) and zero wavelength difference (r/ = 0) as calculated numerically. Dashed lines 
show the corresponding fitting function. Top Left: log(C/) = —0.25, — 0.5, — 0.75, — 1.0 
(from above). Top Right: log{U) = 0.5, 0.25, (from above). Bottom Left: log(C/) = 
0.75, 1.0, 1.25, 1.5 (from above). Bottom Right: log(C/) = 1.75, 2.0, 2.25, 2.5, 2.75, 3.0 
(from above). 
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Fig. 5. — Solid lines show the correlation as a function of wavelength difference t] for zero 
separation (r = 0) and zero source size (r^ = 0) as calculated numerically. Dashed lines 
show the corresponding fitting function. Top Left: log(f/) = —0.25, — 0.5, — 0.75, — 1.0 
(from above). Top Right: log{U) = 0.5, 0.25, (from above). Bottom Left: log{U) = 
0.75, 1.0, 1.25, 1.5 (from above). Bottom Right: log{U) = 1.75, 2.0 (from above). 
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Figure 3 shows that the flux correlation decays with increasing lag. The decay is de- 
scribed by a single lengthscale when U is less than a few (weak scattering regime), whereas 
the decay occurs on two distinct lengthscales for larger values of U (strong scattering). In 
the latter regime, the shorter of the two scales is the diffractive scale and the longer is the 
refractive scale (see Rickett 1990; Narayan 1993; for a review of the underlying physics of 
diffractive and refractive scintillation) , and the ratio of the two scales increases with inreas- 
ing U as C/^/^. The fluctuation power on the diffractive scale remains large for all values 
of [/, whereas that on the refractive scale decays with increasing U. All of these proper- 
ties are well-known and asymptotic results are available in the limit of both very large and 
very small U. The fltting function we have developed handles the asymptotic regimes well, 
but more importantly, it also models the analytically intractable regime corresponding to 
U ~ 1 — 10 quite satisfactorily. The most serious qualitative discrepancy is that the fitting 
function gives a positive value of the correlation for all lags, whereas the numerically com- 
puted correlation sometimes goes negative, especially for [/ ~ 1 (e.g., the top two panels of 
Fig. 3; see also Figs. 6, 7). More pronounced negative "overshoot" is seen in scintillation of 
some extragalactic intra-day variables and has been interpreted as evidence for anisotropic 
scattering(Rickett 2002; Bignall et al. 2003). 

Figure 4 shows the dependence of the flux correlation as a function of the source size r^, 
and Figure 5 shows the dependence of the cross-correlation as a function of the dimensionless 
frequency difference 77. Once again, the results clearly indicate the transition from a single 
scale in weak scattering to two scales in strong scattering. The fltting function does a 
remarkably good job for all values of U and all choices of and 77. 

Figures 6, 7, 8 show the dependence of W{U, r, Vg, 77) as a function of pairs of the three 
variables, r, r^, 77, keeping the third variable flxed at zero. These two-dimensional correlations 
are relatively less well explored in the literature, even in the asymptotic regimes, but our 
results are generally consistent with previous work wherever a comparison is possible. The 
behavior of W{U, r, r^, 0) in Figure 6 is easy to understand. With increasing source size, the 
fluctuations are progressively smoothed, so the fluctuation amplitude decreases and the scale 
of the fluctuations increases, just as one would expect. The dependences in Figures 7 and 8 
are less obvious. With increasing frequency difference r], the decay scale of the correlation 
increases both as a function of r (Fig. 7) and as a function of (Fig. 8); this is especially 
obvious in the plots for weak scintillation and strong diffractive scintillation. Equivalent ly, 
the decorrelation bandwidth increases with increasing r and/or r^. Chashei & Shishov (1976) 
showed such a dependence for strong diffractive scintillation for the particular case of a 
"square law" spectrum, a = 2. No analytical results have been reported for a Kolmogorov 
spectrum a = 5/3, but we see from our numerical work that the overall behavior is similar. 
In any case. Figures 6, 7 and 8 show that the fltting function reproduces the exact numerical 
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results quite well for all values of U and choices of r, r^, 77. 

The behavior of flux scintillation in the r — 77 plane has gained prominence in recent 
years as a result of the discovery of parabohc "arcs" and "arclets" in the so-called secondary 
spectrum (Stinebring et al. 2001; Cordes et al. 2005). Appendix C develops a rudimentary 
analytical theory of arcs that is vahd in the two asymptotic regimes U <^ 1 and U ^ 1. 
The discussion there complements the more detailed physical approach taken by Cordes et 
al. (2005). The behavior of arcs in the transition regime U ~ few, and in the presence of 
anisotropic turbulence, deserves further analytical and numerical study, but it is beyond the 
scope of this work. Researchers interested in arcs should be warned that our fitting functions 
were not designed with arcs in mind. 

We have also compared the fitting function with the numerical results on W{U,r,rs,rj) 
for the case when all three variables, r, r^, 77, are varied, but we do not show the corresponding 
plots. By comparing the fitting function and the numerical results over an extensive grid of 
values of U, r, Vg and 77, we find that the maximum error in this four-dimensional space for a 
gaussian source is 0.047. The maximum error is larger for a tophat source, ~ 0.06, perhaps 
because the simplification of replacing the two source sizes r^i and with a single effective 
size Tg (eq. 18) is less well motivated in that case. 

4.5. Logic Behind the Fitting Function 

Although the fitting function described in the previous subsections was obtained to a 
large extent by a combination of intuition and trial-and-error, we tried to draw on analytical 
clues from Appendix B wherever possible. 

Consider first Wq{U). Equation (B5) shows that in the limit of very weak scattering 
{U <^ 1) the fiux variations have a mean square amphtude equal to 0.7729C/; the first 
term in the fitting function (21) is designed to satisfy this limit. Similarly, for very strong 
scattering ([/ ^ 1), equation (Bll) shows that the refractive fiuctuations have an amphtude 
0.2380C/~°-^ while equation (B14) indicates that the diffractive fiuctuations have amplitude 
1-1-0. 2380f/~° ''. Thus, the total mean square flux variations is equal to 1-1-0. 4760f/^°'^, which 
is ensured by the second term in equation (21) when U ^ 1. Also, we see that the diffractive 
ffuctuations have a baseline amplitude of unity and that the excess ffuctuations above unity 
are divided equally between diffractive and refractive ffuctuations. This is the motivatation 
behind the particular split used in equation (38). The particular functional form that we 
have chosen in equation (21) to interpolate between the weak and strong scattering limits is 
completely arbitrary. In fact, it is easy to find other forms that model the transition region 
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Fig. 6. — Solid lines show the numerically calculated correlation as a function of separation 
r for a series of values of r^, for rj = 0. In each panel, from above, the curves correspond 
to r, = 0, 10-2, 10-l•^ 10-\ 10-°•^ 1, 10° -\ 10\ 10^ '\ 10'\ The dashed lines show the 
fitting function. Top Left: U = 1. Top Right: U = 3.162. Bottom Left: U = 10. Bottom 
Right: U = 31.62. 
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Fig. 7. — Solid lines show the numerically calculated correlation as a function of separation 

r for a series of values of ( (see eq. 27 for the definition), for r,, = 0. In each panel, from 
above, the curves correspond to ( = 0, 10^^, lO^"'^ ''', 10^^, lO^'^ '''^ 1, 10*^'^, 10. The dashed 
lines show the fitting function. Top Left: U = 1. Top Right: U = 3.162. Bottom Left: 
U = 10. Bottom Right: U = 31.62. 
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Fig. 8. — Solid lines show the correlation as a function of source size for a se- 
ries of values of (, for r = 0. In each panel, from above, the curves correspond to 
( = 0, 10"^ 10-^■^ 10-°■^ 1, 10°-^ 10. The dashed lines show the fitting func- 

tion. Top Left: U = 1. Top Right: U = 3.162. Bottom Left: U = 10. Bottom Right: 
U = 31.62. 
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around [/ ~ 1 equally well. 

Consider next the flux correlation as a function of separation r, for a point source 
{vg — 0) and zero wavelength difference (77 = 0). Equation (B13) shows that, for diffrac- 
tive scintillation in the strong scattering regime with a — 5/3, the correlation varies as 
exp(— C/r^/^). When Vg — rj — and C/ » 1, the parameter Ri in equation (24) tends to R5 
in equation (40), whose leading term is proportional to U~^^^. Thus, the r-dependence of 
the fitting function has the correct functional form. Unfortunately, the analytical result (B4) 
does not give a very convenient expression for Fr in the weak scattering limit. However, since 
we require the fitting function to vary smoothly across the transition from strong diffrac- 
tive to weak scattering, we use the same functional form (24) for both regimes. Similarly, 
the expression (B9) does not provide a useful approximation for the factor for strong 
refractive scintillation. By numerical experimentation we determined that F^^r cuts off more 
rapidly compared to F^; we chose the index in (41) to be 7/3 out of a sense of "symmetry." 

Appendix B.3 discusses some asymptotic results for a finite source size. For large source 
size, the flux correlation is shown to decline as a power law in r^, though with a different 
index in the different regimes. In refractive scintillation, the decline is described by oc 
(r^/rref)"^''^ (eq. B24 for a = 5/3) and this is hardwired into the fitting function via 
equation (42). The situation in the case of weak and strong diffractive scintillation is more 
complicated. The former has a scaling W cc Vs (cq. B23) and the latter W oc (rs/rdifr)~^ 
with j3 < 2 (eq. B26). The two indices arc thus different. However, we have required our 
fitting function to have a common set of scalings for weak and diffractive scintillation. (This 
is in the interests of simplicity and for smoothness across U = Uq). Because of this, we 
allowed the index on {r-g/ R2) in equation (25) to be a free parameter and adjusted it to give 
the best overall fit to the numerical data; this resulted in a value for the index of 1.81. 

Finally, consider the correlation as a function of the parameter t] for r = = 0. 
First, notice that r] is defined such that it is equal to AA/A when the two wavelengths are 
close to each other, but it tends to unity as the wavelengths differ by an arbitrarily large 
amount. Because of the latter property, we have found it more convenient to use the modified 
parameter C, defined in equation (27) in the fitting function. This parameter is equal to 77 
when 77 <^ 1, but it is proportional to the wavelength ratio A2/A1 as 77 — > 1. 

In weak scattering and when 77, C, <^1, equation (B5) shows that the correlation varies 
as 1 — 77^/^, i.e., roughly as (1 + C^^^)~^- The functional form chosen for in equation (26) 
is designed to satisfy this hmiting result (note that Oi tends to a constant when U <^ 1). In 

the strong diffractive regime and in the limit of large r/ (large compared to the decorrelation 
bandwidth but still small compared to unity), equation (B20) shows that the correlation 
declines as t;"^?/"^ "^. This asymptotic dependence is ensured by the term a2(t/)C^ in equation 
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(26). Finally, in the strong refractive regime, equation (Bll) shows that the correlation varies 
as 1 — r/^ for small rj. This is taken care of by the term a4(C/)C^ in equation (43), where we 
note that tends to a constant in the limit when U ^ 1. 

While we have tried to respect asymptotic results for [/ <S 1 and C/ » 1 to the extent 
possible, our primary interest is the regime in between where U is neither very small nor 
very large. We have picked functional forms for the various terms in the fitting function 
such that they go smoothly between the two limits and agree as closely as possible with the 
numerical results. This is relatively easy to achieve when only one of the three parameters, 
r, Ts, 1], is varied and the other two arc set to zero (sec the results shown in Figs. 3, 4, 5). 
However, we are also interested in variations of two (e.g., the arc phenomenon. Appendix C) 
or even three of these parameters simultaneously. Finding a reasonable fitting function to 
handle these regimes is more challenging. The particular function described in §4.2 and §4.3 
could doubtless be improved, but it appears to fit the numerical results adequately over the 
entire parameter range of interest. 



5. Summciry 

We provide in §4 of this paper a fitting function for a fairly general correlation, W{U, r, Vg, rj), 
that describes the statistics of flux scintillations of compact radio sources. The function has 
been optimized for two source shapes, gaussian and tophat; the latter may be appropriate 
for gamma-ray burst afterglows (Sari 1998). We also allow for different source sizes and 
observation frequencies for the two flux measurements being correlated, which may again be 
useful for interpreting afterglow observations. 

We have not speciflcally allowed for the flnite integration time or flnite bandwidth 
over which each flux measurement is made. Finite integration time causes any given flux 
correlation to correspond not to a single value of r but to a range of values, i.e., W{U, r, r^, rj) 
is smeared in r by a convolution. Similarly, finite bandwidth leads to smearing in rj. These 
effects can be modeled by convolving the fitting function with the appropriate broadening 
functions in r and rj. Alternatively, they could be incorporated directly into the fitting 
function itself through additional factors (we have not attempted this here). In practice, 
with current technological limits, flnite bandwidth is unlikely to be important in the regimes 
of weak or moderate scattering {U < a, few). 

The results we have presented are for a single thin scattering screen. This is probably a 
reasonable model since the scattering regions in the interstellar medium tends to be clumpy. 
Therefore, the scattering is dominated by one or at most a few distinct screens. The model 
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is also specific to a Kolmogorov spectrum of fluctuations in the scattering medium. This 
again is not unreasonable (e.g., Armstrong et al. 1981), although Boldyrev & Gwinn (2003) 
have recently challenged the usual assumption, which we have adopted, that the statistics 
of these fluctuations are gaussian. 

The most serious simplifying assumption in this work is that we have taken the scattering 
to be isotropic, so that the scintillation correlation function depends on separation r but not 
on the particular direction of r transverse to the line-of-sight. A number of observations 
(e.g., Wilkinson et al. 1994; Dennett-Thorpe & de Bruyn 2003; and references therein) 
have shown that the scattering irregularities in the interstellar medium arc anisotropic, 
presumably because blobs are elongated parallel to the local magnetic field. It would be 
useful to generalize the work described here to include anisotropy. 

We thank Edo Berger, Shri Kulkarni, and Stan Flattc for useful discussions, and our 
anonymous referee for a detailed report containing many very helpful suggestions. This work 
was supported in part by NSF grant AST-0307433. 



A. Derivation of the flux correlation integral 

The goal of this appendix is to derive the expression (A4) for the fourth-order coherence 
function in the special case that the points coincide in two pairs, so that it reduces to a 
flux correlation. Codona et al. (1986) have given a similar result for frequency- independent 
refractive fluctuations, e.g. atmospheric rather than interstellar turbulence. GN89 addressed 
the plasma case but assumed small frequency differences, an approximation that is not made 
here. 



With b = 0, eq. (2.5.4) of GN89 becomes 

1 



r4(0;r;i/i,i/2) 
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exp <|i7r 
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Xiz 
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exp 



"2 ('^(^1' '^l) ~ '^(^3, Al) - 0(X2, A2) + 0(X4, A2)) 



(Al) 



Here X — c/u and z is the effective distance to the thin screen, which impresses a phase shift 
0(x, A) on rays of wavelength A that strike the screen at x. The screen is parallel to the 
observer's plane but separated from it by distance -^screen, while z = (-^g^reen + '^source)"''^' where 
-s^source IS the distaucc from the screen to the source. Since involves the departure of the 
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refr active index from unity due to free electrons all along the path of propagation, we assume 
that the oc A. Also, 0(x, A) is proportional to the column density of electrons, which is 
stochastic. The mean-square phase differences impressed upon parallel paths meeting the 
screen at xi and X2 is 



L>(xi-X2;A) = ([(/.(xi,A)-(/.(x2,A)f> 
L>(xi-X2;A') = (^j^ D(xi-X2;A), 



(A2) 



where (...) is an appropriate statistical average. Let A = (AiA2)^^^, p = (A2/Ai)-'^/^, (pi = 
(pi^i, A), and Dij = D{-Xi — x^. A). Since {4>i4>j) = — ^Dij, where cr^ is the variance of 0, 
the second exponential in eq. (Al) is 



exp [p"'^i3 + + Dr2 + - Du - D23] | . 

Next, multiply (Al) by 



(A3) 



J (ic5 fc-^^Xi 



1. 



Then shift variables x — > x+c, which has no effect on (A3) but multiplies the first exponential 
in (Al) by 

' Xi - X3 X2 - X4 



exp 



27ric- 



Xiz A2-2 

After integration over c, the net effect of these steps is to insert the delta functions 



Xi + X2 + X3 + X4 



Xi — X3 X2 — X4 



Xiz 



X2Z 



into the integrand of (Al). So, only two combinations of the xs are independent. It is 
convenient to choose these as 

q = ^(l + P')(Xl-X3) = i(l+p-2)(x2-X4) 
S = -(Xi +X3 -X2 -X4) , 

because the Presnel {i.e. first) exponential in (Al) then becomes 



exp 



(Al + A2)^ 



q-(s + r) 



After integrating out the two delta functions and expressing the remaining integrals in terms 
of q & s, one has 



in which the structure functions (the Ds) are evaluated at the geometric mean A/A1A2 of the 
two wavelengths, whereas the Presnel scale is defined with the arithmetic mean, eq. (4), and 
rj is the dimensionless frequency difference (5). Equation (A4) agrees with eqs. (17)-(18) of 
Codona et al. (1986), apart from what appear to be minor typographical errors^ and for 
changes to allow for the dependence of the plasma refractive index. 



The bulk of this Appendix is devoted to analytical results for the flux correlation in 
various asymptotic regimes. Since the following presentation is rather dense, we begin with 
an index to the main results: 

• Weak scintillation (C/ < 1): Eqs. (B4), (B5). 

• Strong refractive scintillation: Eqs. (B8)-(B11). 

• Strong diffr active scintillation: 

— Total flux variance (i.e., r = = rj): Eq. (B15) 

— Diffractive correlation at r — but rj ^ 0: Eqs. (B18)-(B19). 

— Diffractive correlation at r Q & rj ^ U^"^^" : Eq. (B21). 

§B.3 discusses the effect of a finite source size. 
^The term 1 + {5k /2k) should read 1 — {5k /2k) and vice versa. 




+ D(s-r?q) + D(s + r7q) -D(s-q) -D(s + q) \, (A4) 




B. Asymptotic Results for the Flux Correlation 
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B.l. Weak scattering (C/ < 1) 

The idea here is to take exp(— G'/2) ^ 1 — G/2 in (8). This is justified because on the 
one hand, \G\ ^ O {Uq") since G^Oass/g^oo; and on the other hand, Uq" <^ 1 unless 
q ^ 1, in which case the rapidly oscillating factor exp(iq- s) suppresses the integral anyway: 

W{ci,rj) ^ (27r)25(q) - ^^/y'c^se^'i-^fls-Tyql^ + ls + Tyql"- |s-q|"- |s + q|"] 

+ 0{U^). (Bl) 

The delta function simply reflects the fact that the mean flux is nonzero. In the remaining 
integral, each term has the form 



oo 



J dse''^-^ |s - a|" = e''^-^ J dte''^-^ |t|" = 2ne''^-^ J e+^Mqt) dt. 







The final integral is divergent, but for —1 < a < —1/2 it would be (Abramowitz & Stegun 
1972, §11.4.16) 

CO 

2n J Jo(gt) dt = -2"+^T\a/2 + 1) sm{na/2)q-''-\ (B2) 



We assume that this formula can in fact be used for all a. Presumably the result could be 
justified by inserting a slowly decreasing smooth function of s into the original integrand 
(8), for example exp(— es^), and then taking the limit e — 0. Such a convergence factor is 
implicit in the Presnel approximation to physical optics since the phase screen is actually 
only a local approximation to a closed surface enveloping the observer and therefore has 
finite area. The final result is 

W{ci; 7], U) ^ 2"+^r^(a/2 + 1) sin(7ra/2) Uq-''-^ [cos{7]q'^) - cos(g^)] 

+0{U'), g^O. (B3) 

The fiux correlation (7) works out to 

W{r;r},U) ^ 1 + 2"r(a/2 + 1)C/ [/(I, r) - /(ry, r)] +0{U^), (B4) 

where f{7],r) = [e'^"/^M(-a/2, 1, irV4r?)] ; 

^ T{a){n\y 
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is a confluent hypergeometric function (Abramowitz & Stegun 1972), and 3? denotes the real 
part. The series for M converges at all z. The "1+" on the first line above represents the 
square of the mean flux; we write AW = W — 1. The correlation at zero spatial separation 
is particularly simple: 

AW(0;r],U) ^ 2°'r(a/2 + 1) cos(7ra/A) {l - rj"'/^) U + 0(U^) (B5) 
^ 0.7729 (1 - r;"/2) U + 0{U^) at a = 5/3. 



B.2. Strong scattering 

B.2.1. Refractive regime 

Even when ?7 ^ 1, the weak-scattering formula (B3) is valid if g <^ ^ref; 

5,ef = U-^l". (B6) 

Actually, one can extend the result to g ~ gref as follows. The integral over s is dominated 
by s ~ where the oscillating exponential cuts it off; in this range, 

G^{1- rj^)Uq'' [a + a{a - 2) cos^ 6] (^)^~" + O {Uq^s''-^) , s > ? , (B7) 

where cos^^ is the angle between q and s. Since we assume a < 2, it follows that G <^ t/g" < 
1, so that we may again evaluate (8) by expanding exp(— G/2) to first order in its argument. 
In contrast to the previous section, however, we retain the prefactor exp(— F/2) because it 
cuts off the refractive spectrum sharply at g ~ gfref^ 

PF,ef(q) ~ 2"+^sin(7ra/2)r2(a/2 + 1)(1 -r/2)t/g2-"exp[-F(g;r;,t/)/2], q < q,,,. (B8) 

Integrating over the direction of q so that exp(q ■ s) — > 27rJo(gr), expressing Jq{z) by its 
power series, and integrating term by term, one finds 

2°-ir(n/2 + l)r(l/a - 1) (1 - //-)[/ 0.9 „w„ 



00 



^ r[(4 + 2n)/^-l] / r , 

^ T{A/a - 1) (n!)2 I ^u^"^ \ ' ^ > 



2 



We have introduced the abbreviation 



Uv = ^(Yir^^i) u- (BIO) 
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Unfortunately the power series on the second hne of (B9) is not any sort of hypergeometric 
function, but it converges for all r if a > 1. In particular, 

AWref (0; V: U) « 0.2380(1 - ri")^ U-^'\ a = 5/3. (Bll) 



B.2.2. Diffractive regime 

If the dimensionless frequency difference rj < C/~^/"^, there is also a diffractive contribu- 
tion to W and W. This comes from q » g'^f ~ s in the integral (8). Consider first 77 = 0. 
Then 

-(F + G) [1 + 0(s/gf-2] . 

Keeping only the leading order term gives 

AVFdiff(r;0,t/) « J ^e" '' J dse'^'^ exp {-Us") . (B12) 

But this is just a Fourier transform followed by its inverse. So to leading order, 

AW^diff(r; 0, U) « exp {-Ur") . (B13) 

We have written AWdiflf rather than W because (B12) applies only to q^ef and therefore 
omits both the refractive contribution (B9) and the square of the mean flux, both of which 
stem from from q < gfref- The total flux correlation is 1 + Al^ref + Al^diff- 

At the current level of approximation (B13), the contribution of diffractive scintillation 
to the flux variance (the correlation at r = 0) is simply unity, i.e., equal to the square of 
the mean flux, whereas the contribution from refraction scales as C/~2(2-°=)/a. We can refine 
the diffractive variance by the following trick. With r = and 77 = 0, the integral (7) is 
symmetrical in s & q: 

-|s + q|" + -|s-q|"-s"-g" 



W{O;0,U) = J dqj e^^^'^ exp 1 1/ 

According to §B.2.1, the contribution from s ^ q is 1 + Al^ref(0; 0, U). Symmetry demands 
that the contribution from s <^ q must be the same: 

A^diff (0; 0,U)^ 1 + Al^ref (0; 0, U) , (B14) 

and therefore the total flux variance in a very narrow frequency band is 



1 = l + 2AWre{{0;0,U). (B15) 
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The factor of two reflects equal contributions from the refractive (s » q) and diffractive 
{s -C q) regimes; but this symmetry should not obscure the fact that the diffractive and 
refractive contributions decorrelate on widely different lengthscales rdifr ~ [/~^/" [eq. (B13)] 
and Tref ~ [eq. (B9)], respectively. 

Let us now consider nonzero frequency differences, rj ^ 0. The most sensitive terms are 
the structure functions C/|s ± Tyqj"; other appearances of rj are smaller by at least 0{rf~°'): 



]-U (|s + r;q|"+ |s-r7q| 



(B16) 



This reduces to (B12) at — 0. With the change of variables 

1/2 / 1 \ 

s = (f) (u + v), (u-v), 

r = (2?7)i/V, and V ^ l{2r])"/^U, (B17) 

the double vector integral can be factored into the product of independent integrations over 
u and V that are complex conjugates: 

APFdiff(r;7;,t/) ^ \j{p.V)\\ (B18) 

oo 

where j{p,V) = ^ ^e"'^"'/'-^"" = j Jo{pu) exp {-^^2 - Vu'') udu. 



[The object /(p, V) is equivalent to the two- frequency, two-position second moment of the 
electric field given by eq. (50) of Lambert & Rickett (1999), who also cite previous work on 
the latter quantity. It is well known that fourth moments of the field, such as W, can usually 
be approximated by squares of second moments in very strong scattering (cf. e.g. Codona 
et al. 1986).] The monochromatic limit 77 — > is recovered by taking F — > while keeping 
(2Vy^°'p — U^/°'r constant. We have no general closed-form expression for f{p,V), but at 
r = 0, 



n=0 



2^/" V- r[2(n + l)/a] ( _^_y-2/oX 
2ar]U'^/^^ n\ \ 2 ) 



(B19) 



n>0 

Looking at the first two terms of this sum, one sees that for small 77, the intensity decorrelates 
as 

W{Q- 0, 0) - W{Q- 77, 0) « 2"r {^^^ cos r^'^'^U 
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in agreement with eq. (21) of Shishov et al. (2003) if one assumes that the scattering, which 
they take to be distributed along the hne of sight, is instead concentrated on a thin screen at 
distance i?/2, and if the term cos should have read cos ^-"^^^"^ therein. The second sum 
in (B19) is asymptotic and shows that 

-2 



Aiydiff(0; V, U) ~ 2-^("-^)/"r2(2/Q; + 1) (ryC/^/") if > C/-^/" 
From either sum, one sees that the decorrelation bandwidth is ~ f/-2/Q 



(B20) 



^2 



Let us now consider the dependence on spatial lag r. For t] < g^^^f, this is approximately 
the same as for = 0, eq. (B13). For 7] ^ g^^^, the term iu^ /2 ^ Vu" in (B18) except 
where the integrand is negligibly small anyway. By omitting iv? /2 from the exponential and 
rescaling the dummy variable u — > one can obtain 

2 



9 



277(t//2)Vc 



(B21) 



where g{z) 



a 



T{2/a) 



Jo{zt)e-^" tdt. 







r{2/a 



r[2(n 

n=0 



l)/a] 



z 

T 



[Mathematically g{z\ is essentially the same as the function Q-p-oij) given in eq. (54) of 
Lambert & Rickett (1999), but the latter is interpreted physically as a pulse- broadening 
profile at zero spatial lag.] To be consistent with the assumption ri ^ f/^^/", one ought 
to use (B20) for AlVdifr(0; f], U) in the above formula, but one might want to use the more 
accurate value |/(0, V)]"^ as given by (B19) instead. 



B.3. Finite Source Size 

Equations (14)-(16) show how a finite source size modifies the flux correlation. For zero 
separation {r — 0), for instance, we have 

W{r = 0, 7^ 0) — / W{(i)qdq, (B22) 
27r Jo 

where refers to the effective source defined in equation (18), and we have made use of 
the fact that the Hankel transform S{q) of a source with a finite size cuts off at g ~ l/rg. 
Using this result, we may estimate the effect of a finite source size in the different scattering 
regimes. In the following, we assume for simplicity that r = 77 = 0. 
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In weak scattering, there is a reduction in the flux variations when the source size is 
greater than the Fresnel scale, i.e., > 1 in our units. Thus, the relevant range of W{q) for 
evaluating the integral in equation (B22) is g < l/r^ < 1. In this regime, W{q) ~ (see 
equation B3), and hence we obtain 

W^weaic(r = 0, r„ ?7 = 0) oc rf-\ r, > 1. (B23) 

In strong refractive scintillation, finite source effects are felt when the source size exceeds 
the refractive scale Tref ~ l/?ref- In this regime again W{q) varies as (see eq. B8), and 
so we obtain 

W^ref (r = 0, r„ 77 = 0) oc (r./r^ef)""', r, > r,ef. (B24) 

Finally consider strong diffractivc scintillation. In the asymptotic limit of very strong 
scattering {U ^ 1), W{q) is constant {W{q) ~ white noise) upto q ~ gaiflf (e.g., Goodman 
& Narayan 1985). For a source size > rdis ~ l/o'dis, this gives 

Wdis{r = 0, r„ 77 = 0) oc (r,/rdiff)-^ r, > r^iS: [/ > 1- (B25) 

However, unless U is extremely large, the next order term in W{q) cannot be neglected, 
and it gives a correction that causes W{q) to decrease with increasing q (e.g., see Fig. 1 of 
Goodman & Narayan 1985, where the decrease is apparent even for U — 10^). As a result, 
the scaling due to a finite source size is modified slightly to 

H^diff(r = 0, r„ r] = Q)(x (r./rdiff)"^, r, > rdiff, (3 < 2, (B26) 

where the exact value of /3 depends on [/; /3 is distinctly less than 2 when U is not very large, 
but it tends to 2 as C/ — > oo. 



C. Pctrabolic eltcs in secondeiry spectra 

The dynamical spectrum of a scintillating source is the flux correlation in a two- 
dimensional plane of time and frequency lag. Its Fourier transform with respect to both 
arguments is called the secondary spectrum S2{fu,ft)- When the interesting range of fre- 
quencies is small compared to the mean frequency, Av <^ u, and when the correspondence 
between spatial and temporal lags is governed by the transverse velocity of line of sight, 
r = v±t, the secondary spectrum is related to the cross spectrum (8) by 

S2{U, ft)-^J'^J^dvW {qtv^ + q'e', 77) e^^^-^ , (CI) 
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in which v±_ = \v±\, and v_l & e' are unit vectors on the sky parallel & perpendicular to 
v±, respectively, while qt = '2nrpft/v± is the component of the wavenumber parallel to the 
motion through the scintillation pattern. The Fresnel scale (4) appears here because W and 
q have been expressed in units where rp = 1: that is, the physical wavenumber is rpQ- The 
integration over q' reduces the 2D cross-spectrum to a ID cross-spectrum. 

Secondary spectra often show a concentration of power along parabolic ridges fi, = ^af^ 
(Stincbring ct al. 2001). Cordes et al. (2005) have given a theoretical explanation of this 
phenomenon, according to which the coefficient a = zX^/{cv\), where z is the effective 
distance of the screen (1). 

The existence of parabolic arcs in strong diffractive scattering can be demonstrated from 
equation (B16), which is equivalent to 



WdisicL; V,U)^ J ds e''^-^ exp 



-it/(|s + r;q|" + |s-r/q|°) 



(C2) 



This is a good approximation to the exact relations (8)- (10) when rfUq°' <^ 1 (even though 
U q°^ may be large) so that the neglected terms in the exponential approximately cancel. For 
large q, the integral over s is dominated by the neighborhood of its branch points. The 
leading-order contribution from s, ri(\ can be obtained by setting |s + ?7q|" |2?7q|" (be- 
cause this part of the integrand remains smooth in that neighborhood) and then expanding 
the remaining exponential to first order in its argument; treating s k, —rjq^ similarly, one has 



Wdiff(q;77,C/) ^ e 



-{U/2)\2m\' 



^ 2(27r)'5(q) + K^^Uq-"-' cos{r]q')e 

where —Ka is the coefficient of g^"^^ in eq. (B2). The delta function must be discarded 
since q is large. When the remainder is used in (CI) to calculate the secondary spectrum, 
the exponential in U\2f]q\°' can be replaced by unity if 

U:$>^u'/"ft, (C3) 

l'V± 

because the complex exponential will then effectively limit rj to small values: 
S2{U,ft)^^K^U—j-{q,+q) + 

(C4)" 

The parabolic arcs can now be recognized in the arguments of these delta functions, which 
are softened but not suppressed by the final integration over the transverse wavenumber: 
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X 




{\M-affy'^'H{\M-af^). (C5) 



Here a = }?zlcv\ as before, and the Heaviside function keeps the argument of the square 
root positive. The Heaviside function cuts off the power abruptly for \ft\ > (|/,y|/a)^/^, 
and the inverse square root concentrates power just inside this cutoff. Both observations 
and numerical simulations confirm this general behavior (Cordes et al. 2005), and this is 
the origin of the arcs. Plausibly the arcs would be more prominent if the electron-density 
spectrum were anisotropic with stronger scattering along than along e', because the 
integration (C4) would have less of a smearing effect. Presumably also, arcs should be more 
prominent when scattering is dominated by a thin "screen" rather than distributed along the 
line of sight, since integration over z would further soften the final square root in eq. (C5). 
These points are made by Cordes et al. (2005). 

In strong scattering {U ^ 1), the conditions for validity of the final result (C5) include 
not only eq. (C3) but also qt = 271 ft/ v± ^ f/^-*^/" = gdifr so that the integration over s is 
dominated hj s <^ q. However, the same result obtains in weak scattering as well provided 
27r/t ^ v±/r-F, since the same cos{riq'^) term appears in eq. (B8) and gives rise to delta 
functions in the 77 — > /i/ transform of the 2D cross spectrum. In fact, regardless of the 
magnitude of the scintillation parameter U, it appears that the arcs are due to weak large- 
angle scattering by inhomogeneities smaller than those that are responsible for the main part 
of the scatter-broadened image (Cordes et al. 2005). 
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